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Abstract 

Charged pattern formation on the surfaces of self-assembled cylindrical micelles formed from 
oppositely charged heterogeneous molecules such as cationic and anionic peptide amphiphiles is 
investigated. The net incompatibility x among different components results in the formation 
of segregated domains, whose growth is inhibited by electrostatics. The transition to striped 
phases proceeds through an intermediate structure governed by fluctuations, followed by states 
with various lamellar orientations, which depend on cylinder radius Rc and x- We analyze the 
specific heat, susceptibility S{q*), domain size A = 2Tr/q* and morphology as a function of Rc 
and X- 
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Many heterogenous molecules, including lipids and amphiphiles, assemble into finite size 
aggregates or micelles 

mm. 

Co-assembly of cationic and anionic heterogeneous 
molecules is a route to create functional biomolecular materials, such as stable vesicles for 
drug delivery and bio-active fibers y, yi. Co-assembled cationic and anionic systems 

are ubiquitous in nature since most biomolecules are heterogenous, charged and strongly as- 
sociating in aqueous solutions. The competition of long-range electrostatic and short-range 
interactions among different co-assembled components may lead to the formation of complex 
aggregates with surface charge heterogeneities. These heterogeneities are expected to play 
a prominent role in the fabrication of functional assemblies and in the self-organization of 
the aggregates. 

The assembly of single component charged amphiphilic molecules is restricted by elec- 
trostatic repulsion 0, [J. For example, electrostatics restricts assembly of charged peptide 
amphiphiles (PA) composed of a hydrophobic block connected to a peptide block that favors 
/3-sheet formation. Single component PAs assemble only when the peptide blocks are neu- 
tral (acidic PAs assemble at high pH and basic PAs at low pH conditions) and/or when the 
salt concentration is excessively high. However, at physiological pH conditions, when the 
molecules are charged, stoichiometric mixtures of acidic (— ) and basic (+) PAs co-assemble 
at one percent concentration into nanofibers of diameter about 6 — 8 nm and a few microns 
length, which form a network that resembles extracellular matrix found in living tissue P|. 
The peptide amphiphile charged end groups are exposed to the surface of the fibers, and 
their amino acids sequences can be designed to promote the growth of bone j^, or neural 
cell differentiation [J]. 

Co-assembled neutral molecules, such as phospholipids and cholesterol in mono- or bi- 
layers, are generally unstable. These neutral assemblies generate domains of different density 
or composition along the surface, which coarsen with time [sj and, eventually, lead to the 
dissolution of the co-assembled structure . Instead, in co-assembled chemically dissimilar 
cationic and anionic amphiphile mixtures, the compositional heterogeneities driven by a 
net incompatibility among the different components are stabilized by electrostatics, which 
favors mixing of the charged components. At low temperature T, the competition of the 
interfacial energy (which favors large domains) and the electrostatic interactions (which 
favors small domains) is expected to result in periodic surface charged domains of size 
A. The cohesive energy \E\ of the domains stabilizes the co- assembled periodic structure 



2 



at low T. In analogy with ionic crystals of ions of valence Z and periodicity A at T = 
where E ~ —Z'^/A since in our system Z ~ A^, the energy per particle ~ A^ 
scales as E/N ~ —A leading to robust periodic structures. In this letter we study the 
formation of charged domains in cylindrical aggregates with Coulomb interactions at various 
degrees of incompatibility and cylinder radii. We analyze charge correlations and domain 
orientations, as well as the thermodynamics of the local segregation process. We show that 
phase segregation on the surface of co-assembled cationic-anionic micelles can be readily 
observed when the net incompatibility among dissimilar molecules generated by chemical 
and/or molecular sizes differences is larger than a couple fc^T's. 

We consider each aggregate as a stable structure, thus constraining the cylindrical ge- 
ometry of the micelles and examining only equilibrium surface structures. A binary elec- 
troneutral charged lattice fluid consisting of positive and negative units of equal absolute 
charge \Q+\ = \Q-\, confined to a cylindrical monolayer of radius Rc and length Lc = L^, is 
placed in the center of the box L x L x along the z— axis. The surface of the cylinder 
is filled with spherical units of diameter cr = 1, in such a way, that all units are placed into 
the knots of a triangular lattice of period a = a (Fig. QJ. The net degree of compatibility in 

ksT units (x = n[e+^ — {e + £:++)/2], where n = 6 is the number of nearest neighbors) 

in our incompressible model is simply x = ~^^++, because for simplicity we only consider 
short-range attraction among positively charged molecules (in a compressible model all the 
interactions have to be included). It should be noted here that x describes the effective 
interaction among the co-assembled macromolecules and not among single monomers. In 
the systems with only van der Waals interactions x is proportional to 1/T, however, sys- 
tems with hydrophobic and hydrogen bonding interactions, such as in co-assembled PA's, 
the pair-interaction may have a complex T dependence. Excluded volume and electrostatic 
interactions are also considered. The Hamiltonian of the system reads as follows. 



where Qi is a charge of the i—th unit, eo is the average dielectric constant of the media in 
units e'^/Ana, Vij = |r^- — rj|/cr is a dimensionless distance between i—th and j— th units, rj 
determines the position of the i—th unit in the space, = 1 and the summation in the 
second term is only over nearest neighbors. To exclude interaction with displaced image- 
cylinders, periodical boundary conditions are applied only in the 2;— direction and, thus, 




(1) 
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Lekner summation technique is used to calculate the electrostatic energy of the single 
cylinder system 

Eel _ 1 v;^ v;^ QiQj 



where pij = [xf^ + yfj\ ^^"^ and for m = the terms with i = j are omitted. We report 
standard canonical Monte Carlo simulations following the Metropolis scheme for various 
values of G [—4.5,-0.5], Cq = 10, Rc/cr G [1,8] and Lc/cr = 100. Simple moves in 
the phase space are performed by exchange of two randomly chosen particles. Each system 
is equilibrated during 10^ MC steps per particle and another 10^ MC steps are used to 
perform measurements. The equilibration process is accompanied by a gradual decrease of 
temperature (temperature annealing) from T^ax = 10 to Tmm = 1. To analyze equilibrium 
properties we calculate the heat capacity Cy and the static structure factor 

2, 



N+ 



i,j=0 
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where Sj is a two-dimensional cartesian surface vector and S{q) is averaged among differ- 
ent directions. We focus our attention on the size of the segregated domains and on the 
susceptibility, a degree of correlations in the system. 

The transition from the isotropic to the striped phase begins with the appearance of 
segregated domains (Fig. Q (a)) and proceeds through an intermediate locally correlated 
state governed by fluctuations (Fig. ^(b)). In the planar (2D-monolayer) case, in the limit 
of zero charge, a macroscopic phase segregation in the mean field approximation is expected 
at a critical point x'^^ — 2, as predicted by a simple coarse-grained free energy functional of 
the neutral system, (l^ 

AF(0) /• /(0) + «:|V0|2 



r1 



-d\ (4) 



where 



is the local charge concentration, / {(f)) = (f)\n(f) + (1 — (j)) \n (1 — (f)) — , 
K ^ jl4|. In the charged case there is no critical point. Instead, a broad peak 

in the heat capacity (Fig. at a characteristic value denoted by x^^^ is observed. The 
charges restricts the possibility of macroscopic segregation and induces instead favorable 
fluctuations of wave length A. If the charge fluctuations are small, a density variation till 
second order in the free energy in terms of the Fourier components of the density (pq gives 



AF/kBT ~ J2q4'q(l>-q/i'^So{q)), where 5'o"^(g) = 4 - 2% + Xq'^ + Svr/geo. The competition 
between electrostatic, Svr/geo, and gradient, xi'^y energies results in the formation of most 
favorable fluctuations of finite wave length A = 27r/g* and the appearance of a peak in 5*0(5) 
at q* = (Att / toXY^'^ shown in the computed S{q) in Fig. 01 (a). This scaling regime is 
possible for x < X^^^ (fluctuations are expected to modify this mean field scaling) (Fig. El (c)). 
On the other hand, in the high x limit (low T), strongly segregated charged domains of size 
A with a well defined line tension 7 ~ x develop. The competition of the interfacial 
energy among the segregated domains, 7A, and the electrostatic penalty associated with 
creating a charge domain, (A^)^/(eoA), gives another scaling limit A ~ i^oXY^"^ Ql- The 
second scaling regime is found to exist for x > X^'^\ where x^"^^ is defined below. 

The mean-field analysis of Eq. HI erroneously predicts a continuous transition to a periodic 
structure signaled by a mean field susceptibility So{q*) = cxd at Xcp = 2 + (3/2) [An y/Xcp/^o] "^^^ ■ 
It is well known that the transitions from isotropic to periodic structures cannot be con- 
tinuous in 3D In 3D the electrostatically driven microphase separation ^| is often 



suppressed by various effects [l8|. In 2D, since there is no long range order jl9|, the possibil- 
ity of a classical thermodynamic transition to periodic structures is questionable In ID, 
on the other hand, fluctuations destroy the possibility of either macrophase or nano-phase 
transitions at non zero T [2^. We expect to recover the 2D limit only when Rc 00 and 
the ID case when Rc ^ 0. 

Appearance of large domains on the surface of the cylinder of final radius applies geo- 
metrical restrictions on the size and symmetry of the domains, especially in the ID-limit. 
These geometrical restrictions may influence physical properties, such as the susceptibility. 
Figure El (b) shows a non-monotonic dependence of S{q*) with x- We find that S{q*) is 
independent of Rc for x < X^'^\ while for x > X*^^^ it shows a non monotonic dependence 
on Rc. The point x^"^^ corresponds to the cross-over transition from the state, where the 
segregated domains are locally correlated, to the striped state, where the domain size scales 
as A ~ (coX)^^^ s-nd the domains morphology depends on R^. The melting of striped struc- 
tures in the limit Rc ^ 00 may occur via the Kosterlitz-Thouless mechanism 2M- Besides 
the appearance of dislocations and other defects (Fig. ^(c)), the striped state on the cylin- 
der is richer due to the possibility of different symmetries. We find a broad peak in S{q*) 
versus Rc (Fig. 131(a)) for each x > X^'^\ implicitly denoted by x^^\ that moves to larger Rc 
values with increasing x- For intermediate x > X^"^^ ^he structures are strongly dependent 
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on Rc. For very wide cylinders {q*Rc > 2.5) we find a defect mediated striped state. Witfi 
decreasing Rc the striped state turns into tlie spiral/zigzag state (1 < q*Rc < 2.5). Witli 
furtfier decreasing R^ the spiral/zigzag pliase, first, turns into a parallel lamella and then 
into the ring phase {q*Rc < 1). With increasing the energy x, the spiral state turns into the 
zigzag state and then into the parallel striped phase. Our simulation shows that for large x 
the parallel stripe phase dominates unless the q*Rc ~ 1. 

We calculate an order parameter with the purpose of characterizing the geometry of the 
domains and to find their preferable morphology. We begin by defining independent clusters 
on the surface of the cylinder. For each cluster C we calculate the inertia matrix with 
components 

. Nc 

where Nc is the number of ions in the C— cluster and r° is the a— th cartesian component 
of the position vector of the i—th ion. To characterize the aniso trop y of the cluster we use 
the "asphericity parameter" introduced by Rudnick and Caspar 23 1 

^ /(rr|r|)--3M|r|\ 

\ iTrlT]f I; 

where (i?^, i?|) are the eigenvalues of the matrix Tap which defining the three principal 
axes of the cluster, Tr[T] = Rj + R^ + R^ is the trace, M[T] = RjRj + RjR^ + R^Rl is 
the sum of three minors and A is averaged among all clusters. For spherically symmetric 
objects (-Ri = R2 = R3) A is equal to zero, for thin toroidal objects = R2, R3 <^ -Ri) it 
is equal to A ~ 1/4 (in our case, ring phase) and for long thin stripes (-R2 ^ Ri, R3 ^ Ri) 
y4 ^ 1 (parallel lamella phase). Figure ID (b) shows A as a function of x- For small values 
of X; when system is isotropic, the asphericity parameter takes value A ~ 0.6 independent 
on Rc- For wide cylinders {q*Rc > 1) A grows with x aiid approaches unity at large x 
values. For relatively thin cylinders, q*Rc — 1.5, it approaches unity faster than for cylinders 
with q*Rc — 2.5, due to the fact that A does not distinguish stripes in parallel, spiral or zigzag 
states. Meanwhile, for the stripes disturbed by the defects {q*Rc > 2.5) the asphericity 
parameter takes smaller values than unity. In the case of thin cylinders {q*Rc < 1) the 
asphericity parameter A indicates formation of ring stripes. The transition from the parallel 
lamellar phase to the ring phase is found at x denoted by x^^^ ■ 

Our results are summarized in Fig. HI (c), where we show a schematic {x,Rc) phase 



diagram. We predict various nano-phases obtained as a result of the interplay of short- 
and long-range forces among stoichiometric charged components confined on the surface of 
a cylindrical fiber. Our study suggests that besides modifying the size and the chemistry of 
the molecules to induce different domain symmetries, one can co-assemble basic and acidic 
molecules with different total charge per molecule, e. g. + 2 and —1, which may induce other 
domain symmetries along the surface of the cylinder jlO| . Various different symmetries and 
interesting phenomena appear as a result of the interplay between the symmetry of the 
surface domains and the geometry of the aggregates. The constraint of tiling the surface 
of aggregates with finite size domains of different symmetries is a topic of interest and is 
relevant to studies of viral symmetry and other self-assembled supra-molecular structures. 

We acknowledge valuable discussions with S. I. Stupp, F. J. Solis and A. Kudlay. This 
work is supported by NSF grant numbers DMR-0414446 and DMR-0076097, and by NIH 
grant number GM62109. 
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FIG. 1: Snapshots of typical configurations of (a) isotropic phase (x = 4.5), (b) locahy correlated 
state (x = 6.75), (c) striped state with defects (x = 9.0) and (d) parallel striped phase (x = 12.0). 
Rc/a = 8. 
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FIG. 2: Dimensionless, normalized by A^, heat capacity as a function of the pair attraction energy x 
for different Re- 
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FIG. 3: (a) Static structure factor for Rc = 5 for different x, (b) magnitude of the peak of the 
structure factor, S{q*), as a function of x for different Rc values, (c) double logarithmic plot of the 
peak position q* versus x for different Rc- 
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FIG. 4: a) Susceptibility S{q*) as a function of Rc near extremum, b) asphericity order parameter 
as a function of x and Rc, c) schematic phase diagram (x, Rc) implied by results of simulation, 
where I - isotropic phase, II - parallel lamella phase, _L— ring phase, LC - locally correlated state, 
S - helically twisted state, Z - zigzag state, DM - defect mediated state. Lines (1), (2), (3) and 
(4) correspond to x^^\ X^'^\ X^^^ X^^'^ ■ 
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